arXiv:1502.06884v2 [gr-qc] 5Jun2015 


Black hole initial data without elliptic 

equations 

Istvan Racz and Jeffrey Winicour ^ 


§ Wigner RCP, 

H-1121 Budapest, Konkoly Thege Miklos ut 29-33., Hungary 

b Department of Physics and Astronomy, 

University of Pittsburg, Pittsburgh, PA, 15260, USA 

Abstract 

We explore whether a new method to solve the constraints of Einstein’s 
equations, which does not involve elliptic equations, can be applied to provide 
initial data for black holes. We show that this method can be successfully 
applied to a nonlinear perturbation of a Schwarzschild black hole by establishing 
the well-posedness of the resulting constraint problem. We discuss its possible 
generalization to the boosted, spinning multiple black hole problem. 


1 Introduction 

The prescription of physically realistic initial data for black holes is a crucial ingredient 
to the simulation of the inspiral and merger of binary black holes and the computation 
or the radiated gravitational waveform. Initialization of the simulation is a challenging 
problem due to the nonlinear constraint equations that the data must satisfy. The 
traditional solution expresses the constraints in the form of elliptic equations. Here 
we consider a radically new method of solving the constraints which does not require 
elliptic solvers [T]. We show, at least for nonlinear perturbations of Schwarzschild 
black hole data, that the Hamiltonian and momentum constraints lead to a well-posed 
strongly hyperbolic problem whose solutions satisfy the full constraint system. The 
possibility of extending this approach to binary black holes offers a simple alternative 
way to provide boundary conditions for the initialization problem that might prove 
to be more physically realistic. 

The inspiral and merger of a binary black hole is expected to be the strongest 
possible source of gravitational radiation for the emerging field of gravitational wave 
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astronomy. The details of the gravitational waveform supplied by numerical sim¬ 
ulation is a key tool to enhance detection of the gravitational signal and interpret 
its scientific content. It is thus important that the initial data does not introduce 
spurious effects into the waveform. Such “junk radiation” is common to all current 
methods for supplying initial data and appears early in the simulation as a high fre¬ 
quency component of the waveform. This can be a troublesome feature with regard 
to matching the waveform in the nonlinear regime spanned by the simulation to the 
post-Newtonian chirp waveform provided by perturbation theory. The initial param¬ 
eters governing the black hole spins, mass ratio and ellipticity of the binary orbit have 
to be adjusted to include the effect of this transitory period. As a result, it becomes 
difficult to match exactly to the parameters governing the post-Newtonian orbit. In 
addition, although the high frequency component of the junk radiation appears to 
dissipate after some early transitory period, there is no quantitative measure of its 
low frequency component which might affect the ensuing waveform. 

All initialization methods presently in use reduce the constraint problem to a 
system of elliptic equations, which require boundary conditions at inner boundaries 
in the strong field region surrounding the singularities inside the black holes, as well 
as at an outer boundary surrounding the system. The new method we consider here 
only requires data on the outer boundary, which is in the weak field region where 
the choice of boundary data can be guided by asymptotic flatness. The constraints 
are then satisfied by an inward “evolution” of the hyperbolic system along radial 
streamlines. 

The initial data for solving Einstein’s equations consist of a pair of symmetric 
tensor fields ( hij,Kij ) on a smooth three-dimensional manifold E, where is a 
Riemannian metric and K t j is interpreted as the extrinsic curvature of E after its 
embedding in a 4-dimensional space-time. The constraints on a vacuum solution (see 
e.g. Refs. PIS]) consist of 

(3) R + (K j j ) 2 - KjjK ij = 0 , (1.1) 

DjK j i — DiK j j — 0 , ( 1 . 2 ) 

( 3 ) 

where R and Di denote the scalar curvature and the covariant derivative operator 
associated with hjj , respectively. 

The standard approach to solving the constraints is based upon the conformal 
method, introduced by Lichnerowicz [4] to recast the Hamiltonian constraint (HU 
as an elliptic equation and later extended by York ®E] to reduce the momentum 
constraint (11.21) also to an elliptic system. For a review of the historic implementation 
of this method in numerical relativity see [7]. 

A major obstacle in prescribing black hole initial data is the presence of a sin¬ 
gularity inside the black hole. The initial strategy for handling the singularity was 
the excision of the singular region inside the black hole [8]. In this case, an artificial 


2 




inner boundary condition for the elliptic system is posed on boundaries inside the 
apparent horizons surrounding the individual black holes. Other strategies have since 
been proposed. One is the puncture method in which the initial hypersurface extends 
though a wormhole to an internal asymptotically flat spatial infinity, which is then 
treated by conformal compactihcation |Qj. Here the freedom in the choice of confor¬ 
mal factor governing the compactihcation enters as an effective boundary condition. 
In addition, it is known that the puncture quickly changes its nature. In fact, early 
attempts to simulate binary black holes failed until it was realized that the punctures 
must be allowed to move. Studies of this feature in the case of a single black hole 
revealed that the puncture quickly transits from the internal spatial infinity to an 
internal timelke infinity [lOj. This realization has given rise to the trumpet version 
of initial data, in which the initial Cauchy hypersurface extends to an internal time¬ 
like infinity with asymptotically finite surface area maim. Trumpet data offers a 
promising alternative to puncture data but its merits have not yet been extensively 
explored in binary black hole simulations m 

Coupled to these techniques for avoiding singularities is the choice of initial time 
slice. For example, there are many ways to prescribe Schwarzschild initial data de¬ 
pending, say, upon whether the initial Cauchy hypersurface is time symmetric or 
horizon penetrating. Here we will focus on initial data in Kerr-Schild form mm, 
which for the Schwarzschild case corresponds to ingoing Eddington-Finklestein coordi¬ 
nates, which extend from spatial infinity to the singularity and penetrate the horizon. 
The new approach to solving the constraints that we consider becomes degenerate 
for a time symmetric initial slice, whose extrinsic curvature vanishes. However, time 
symmetric space times contain as much ingoing as outgoing gravitational waves, so 
they are not the appropriate physical models for studying binary waveforms. Al¬ 
though our focus here is on data in Kerr-Schild form, we do not wish to imply that 
this approach would not work for puncture or trumpet data. 

A very attractive feature of Kerr-Schild initial data is that it provides a preferred 
Minkowski background to construct boosted black holes by means of a Lorentz trans¬ 
formation. Two independent ways of prescribing Kerr-Schild initial data have been 
proposed. In one version, the 4-dimensional aspect of the Kerr-Schild ansatz is pre¬ 
served as much as possible [15]. This leads to a workable scheme for superimposing 
non-spinning black holes but the generalization to the spinning case remains prob¬ 
lematic. In the other case, the Kerr-Schild ansatz is loosened to a 3-dimensional 
version that allows superposition of multiple spinning black holes [16]. This has been 
implemented to provide data for boosted, spinning binary black holes and plays an 
important role in current simulations H3- 

There are several variants to the new method of solving the constraints proposed 
in HJEHUIM, depending upon which components of the initial data are assigned 
freely. They all avoid elliptic equations. Here we apply the simplest of these variants 
to the initial data problem for black holes. In this variant, the Hamiltonian and 
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momentum constraints constitute a strongly hyperbolic system which only requires 
data on a 2-surface surrounding the black holes. 

In Sec. El we review this new approach. In Sec. [3] we show that the require¬ 
ments for well-posedness of the underlying algebraic-hyperbolic constraint problem 
are satisfied by a Schwarzschild black hole described in Kerr-Schild form. In Sec. El 
we present an explicit proof that nonlinear perturbations of Schwarzschild black hole 
data in Kerr-Schild form lead to a well-posed strongly hyperbolic problem. 

In Sec. El we conclude with a discussion of the possibility of extending this ap¬ 
proach to general data for a system of boosted, spinning multiple black holes. We 
show how the initial metric data for multiple black holes can be freely prescribed in 4- 
dimensional superimposed Kerr-Schild form for the individual boosted, spinning black 
holes. Two pieces of extrinsic curvature data, which represents the two gravitational 
degrees of freedom, can also be freely prescribed by superimposing the individual 
black hole data. The remaining extrinsic curvature data is then determined by the 
algebraic-hyperbolic constraint system. In a linear theory, the superposition of such 
non-radiative data would lead to a non-radiative solution. This suggests that this new 
method may offer an alternative approach to suppressing junk radiation and to con¬ 
trolling the effect of initial data on a binary orbit. However, due to the nonlinearity 
of Einstein’s equations, there is no guarantee that, in the strong field region between 
the individual black holes, this superimposed free data does not introduce spurious 
radiation. A completely analytic resolution of these issues does not seem possible. 
A major motivation for this paper is to encourage the numerical experimentation 
necessary to explore the merit and feasibility of this new approach. 


2 A new approach to the constraints 

We assume that the topology of £ allows a smooth foliation by a one-parameter family 
of homologous two-surfaces. In the application to black hole initial data, we assume 
for simplicity a foliation 5? p by topological spheres described by the level surfaces 
p = const of a smooth function. 

Choose now a vector field p l on £ such that p l dip = 1. Then the unit normal h l 
to Sf p has the decomposition 


n* = N 1 [p i -N i ], (2.1) 

where the ‘lapse’ N and ‘shift’ N l of the vector field p l are determined by fq = Ndip 
and N l = 7 * j p 7 , with 7 1 j = 8 l j — h l hj. 

The 3-metric hq on £ then has the 2 + 1 decomposition 

h{j 7 ij T 8ipij , (2.2) 
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where 7 ^ is the metric induced on the surfaces 5? p . The extrinsic curvature K V} of 
Sf p is given by 

Kij = 7 \ Di hj = \ &hlij ■ (2.3) 

The extrinsic curvature K tJ of E, which forms part of the initial data, has the 
decomposition 

= k hihj + [hi k j + hj kj] + K L j , (2-4) 

where k, = h k h l Km, kj = 7 h l K ki and K^- = J k i A /j K k i_. Here we use boldfaced 
symbols to indicate tensor fields tangent to 5^ p . In addition, we shall denote the 
trace and trace free parts of K l3 and K, :j by K\ = 7 kl Km, K l i = 7 fci K K tJ = 

Kij - \ 7 ij K\ and K tj = K tJ - \ 7 tj K\, respectively. 

By replacing the initial data set ( hij , K i3 ) by the seven fields (IV, N\ 'jij, K ?J -, k, kj, K(), 
the Hamiltonian and momentum constraints (11.11) and (11.21) can be expressed as pQ 
(see also dHUBEQ]) 

^(K \) - D l kj + 2 h l kj - [ k - i (K'j) ] {K\) + K u K kl = 0 , (2.5) 

2^kj + (K^-^k A(K Z j) - 2k z Aki] + (2K l i )- 1 A f ] R - K kl K kl ] 

+(k\) kj + [« - i (K',) ] hi - h l K/j + D'Kjj = 0 , (2.6) 

where k, is determined by 

#c = (2 K ',)" 1 [ K fcI K w + 2 k'kj - | (K'j ) 2 - < 3 ) i? ], (2.7) 

-Dj and R denote the covariant derivative operator and scalar curvature associated 
with 7 jj, respectively, and h *. = h l D[hk = — Dk(lnN). Here (12.71) provides an al¬ 
gebraic solution to the Hamiltonian constraint ( 11 . 11 ) (for more details see pQ). The 
four quantities (n, kj, K \) are subject to the constraints whereas the remaining eight 
variables (N, N l , ^ij, Kjj) are freely-specifiable throughout E. Here K ^ encodes the 
two free gravitational degrees of freedom. 

Given the free data (N, N\ 7 ^, K,j), the equations (I2.6D - (I2.7D were shown to com¬ 
prise a first order strongly hyperbolic system for the vector valued variable (K^,kj) 
provided n and K \ are of opposite sign, 

kK\ = -C 2 , C ^ 0. (2.8) 

It was also verified in [I] that, given the values of (kj, K \) on some “initial” surface J ^ 0 
satisfying (12.81) . solutions to the nonlinear system (I2.5D - (I2.7D exist (at least locally) 
in a neighborhood of 2 /q, and that the fields ( hij,Kij ) built up from these solutions 
satisfy the full constraint system (ll.ll) - (ll.2D . 
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3 The free and constrained Schwarzschild data 


The successful application of this new approach to the constraint problem depends 
upon a judicious choice of gauge, determined by the lapse of the initial Cauchy hy¬ 
persurface E, and a judicious choice of foliation 5^ p . We begin by considering data in 
Kerr-Schild form, in which the space-time metric has the form 

g ab = Vab + 2 HUb , 9 ab = V ab - 2 Hi a t , (3.9) 

where H is a smooth function (except at singularities) on M 4 and i a is null with 
respect to both g ab and an implicit background Minkowski metric rj ab . In inertial 
coordinates (f, x l ) adapted to r] a b, 

g ab dx a dx h = (-1 + 2 H£ t 2 )dt 2 + Am t £ i dtdx i + (d lJ + 2 H£ i £ j )dx i dx j , (3.10) 

where £“ = g ab £ b = g ab £ b and g ab £ a £b = V ab £ah = — (£t) 2 + = 0. The Kerr-Schild 

metrics also satisfy the background geodesic condition 

v b %d b £ a = 0 (3.11) 

and wave equation 

if b d a d b H = 0. (3.12) 

We can relate the Kerr-Schild metric to the 3 + 1 decomposition of the space-time 
metric 

Qab h ab n a n b , (3.13) 

where n a is the future directed unit normal to the t = const hypersurfaces. Choose a 
time evolution field t a satisfying t a d a t = 1. Then n a has the decomposition 

n a = N-\t a - N a ), (3.14) 

where N and N a denote the spacetime lapse and shift, determined by 

N — — (t e n e ) ,n a — —Nd a t and N a = h a e t e , (3.15) 

respectively. 

In the Kerr-Schild spacetime coordinates (f,x*), the metric has components 


9a0 


—N 2 + NiN* Ni \ 

Nj hij J 


(3.16) 


It follows that 


hij = Sij + 2 Hlilj , h ij 


6 ij - 


2 Htl j 
1 + 2 Hi 2 ’ 


(3.17) 
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(3.18) 


N = 


y/1 + 2 HI 2 ’ 

Ni = 2Hl t li , IV* = 2HN 2 £ t £ i . 

A direct calculation of the extrinsic curvature 

K t j = |ifAi = (27V)- 1 [dAj - (A A + A-A)] 

gives 

N-'Kij = -^ [diiHij) + djiHit)] + N~ 2 d t {Hi i £ j 
+ 2H£ t £ k d k (H£ i £ j ) - H{£{d 3 t t + ^ t ). 


(3.19) 


(3.20) 


For a Kerr spacetime 


H = 


rM 


(3.21) 


(3.22) 


r 2 + a 2 cos 2 6 ’ 

where the Boyer-Lindquist radial coordinate r is related to the Cartesian inertial 
spatial coordinates x l = (xi,X 2 ,x$) according to 


r 2 = 


(p 2 — a 2 ) + \J {p 2 — a 2 ) 2 + 4a 2 x| 


with 

and 


4=1 


2 2 , 2,2 

p = X x + X 2 + x 3 


rxi + ax2 rx2 — aaq £3 


y + a 2 


r 2 + a 2 


(3.23) 


(3.24) 


(3.25) 


As H and £ a are t -independent and £ t — 1, the extrinsic curvature (13.211) simplihes 


to 


/A = -4A^ [4 {Htj) + <9,-(#4) + 2H£ i £ j £ k d k H] . 

Formally, for the purpose of applying the approach in Sec. [2] to a generic inspiral 
and merger, it would be sufficient to show that the required sign condition (12.81) holds 
for a boosted Kerr black hole. Here we restrict our investigation to the Schwarzschild 
case, where the choice of foliation is guided by spherical symmetry and the alge¬ 
braic simplicity allows a clear exposition of the approach. 

For a Schwarzschild black hole, the spin parameter a = 0 and the Kerr-Schild 
form of the metric simplihes to 


H — — , £%= '— = d t r , r 2 = 5^x\Xj , 


(3.26) 
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with lapse 


N = (1 + 2 H)~ 1/2 


(3.27) 


and 3-metric 

hij = 8ij + 2 Hlilj . (3.28) 

(Here — i a is a future directed ingoing null vector, which corresponds to the convention 
for ingoing Eddington-Finklestein coordinates.) Thus 

diH = Xi , dj(H£i) = ^ \r 2 Sij - 2x i x j \ (3.29) 


and (13.261) reduces to 


Ka = 


2 M 


y/l + 2 H 


(S ij -[2 + H]£ i e j ) • 


(3.30) 


We choose the foliation 5^ p by setting p — r, with p l = t, corresponding to the 
; ‘spatial” lapse and shift 

N =Vl + 2H, N* = 0, (3.31) 


unit normal 


•hi = Vl + 2 H li , h l = h VJ hj = 
and intrinsic 2-metric 


/J y/1 + 2 H 


f, 


(3.32) 


7 ij = h ij - nfij = $ij - tdj > 7 4J = ^ ^ • (3-33) 

A straightforward calculation gives the extrinsic curvature components of £, 


2M (1 + H) 

hl = n n K kl = , 

r 2 (1 + 2 Hf' 2 

(3.34) 

ki = 7 k ih l K ki = 0 , 

(3.35) 

Afc A , _ 2M 

W 7.7j^i r 2 yi^ 2 i 7 7 ^’ 

(3.36) 

1 r 2 Vl + 2H 

(3.37) 

w- 

W 

c*. 

'■o. 

I 

toll— 1 

^o. 

o 

(3.38) 


and 


Note that hl and K \ are globally non-vanishing and have opposite sign, in agreement 
with the condition (I2.8j) for strong hyperbolicity. 
















From (12.3p along with 


^hlij = h k d k *{ i:j + 7 kj{dih k ) + 7ik{djfi k ) 

1 


k 

x a 
- O k 


\/l + 2 H r 

XiXfc 

Oik -V 


OC 2 OC j 


X k X j\ rj 

1 

i— 1 

r 2 ) a > 

Vl + 2 H r _ 


d-i 


X 


Vl + 2 H r 


r x/1 + 


2 M 


(4, - *2) 


the extrinsic curvature of the p = r = const foliated surfaces is given by 

1 


K-ij __ ^— | n TT r Yij ; 


VI+ 2 H 


so it follows that 


and 


K\ = 7 kl K kl = 


r ^/1 + 2 H 

Kij = Kij - \ t/ij k\ = 0 . 


(3.39) 


(3.40) 

(3.41) 

(3.42) 


4 Nonlinear perturbations of a Schwarzschild black 
hole 

Here we investigate nonlinear perturbations of the Kerr-Schild initial data for a 
Schwarzschild black hole. In doing so, we simplify the discussion by assigning Schwarzschild 
values to the freely specifiable variables (N,N\ 7 ^, Ky). As a result, the initial 3- 

o /s ^2) 

metric h l3 retains its Schwarzschild value and, in particular. K ?J = 0 and N and R 

* * . O 

have no angular dependence. For a more general perturbation, (. N, N l , 7 ^, Ky) would 
enter as explicit terms in the resulting system for (k, K ; /, k ? ,). 

In this setting, fl2.5l) - fl2.6D reduce to 

=§4(K l i) - D% * [ K - \ (K\) ] (k\) = 0 , (4.43) 

* + (K^-^k A(K'/) -2 k'AV] + (/V) k, = 0, (4.44) 

where k, determined by ( 12 .7j) , reduces to 

« = (2K^)- 1 [2k'k l - 1 (K'O 2 - ( 3 , i?] . (4.45) 

It is easy to check that these equations hold for a Schwarzschild solution, for which 
{3) R = r 4 (i+oj /)2 1 = yi+oj/ V k 2 ; = 0 and neither K z z nor k. have angular depen¬ 

dence. 
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(4.46) 


In spherical coordinates x 1 = (r, x A ), x A = (#,</>), 

7 ijdx l dx J = r 2 qABdx A dx B 

where Qab is the unit sphere metric. Then (14.431) - (14.441) become 

-=1=9, K', - D B k B -[k-\ (K 1 ,) ] (A",) = 0, (4.47) 

1 = drk A + (K'O-'lKSjtK 1 ,) - 2 k B D^k B ] + (A 1 ,) k B = 0 . (4.48) 

y/1 + 2 H 

Now consider nonlinear perturbations of Schwarzschild. We denote by 5V = 
V — Vs the deviation of a variable V from its Schwarzschild value Vs- Then (I4.47D - 
(I4.48D take the form 

1 a BC 

— 7 ====a r aK‘, - -^dcSks = U, (4.49) 

1 /c 2a BD \<n 

dr dk A + —- d A^B = F A , (4.50) 

vl + 2 H K 7 


where F\ and Fa represent lower differential order terms. This is a coupled quasilinear 
system for the vector valued variable U a = (wi, ua) = (5K l i, clk^). The system il 1. 1! Hi ¬ 
ll l.'iOD has matrix form 

d T U a = Cj c d c Up + F a , (4.51) 

where d T = (1 + 2H)~ 1 Fd r , F a = (F { , Fa) and 


£x 1C = 0 , Cr=^q 


BC 


BC 


C ic 


K, 


A ~ 


iS c c 

o A , X* 


K l A 


BC 

A ~ 


r 2 K l , 


n S' 


D u A ' 


(4.52) 

(4.53) 


The requirement that (14.511) is a strongly hyperbolic system |2TI. 22] is that there 
exists a positive bilinear form such that C(u)p a = Ft^C ^ 0 loq is symmetric for 
each choice of uq- It is straightforward to check that such a symmetrizer is given by 


H ii 

H A i 



Hia = 0 , 


H A b = r 2 q A B- 


The positivity of the symmetrizer for perturbations of Schwarzschild, 
H a pv a v^ = - —(u 1 ) 2 + — k A v 1 v A + r 2 q A BV A v B > 0, v a ^ 0, 

K K 


(4.54) 

(4.55) 


(4.56) 


10 































follows from the near Schwarzschild approximations 


K l t _ 2(1 + 2 H) 
k, 1 + H 



( 4 . 57 ) 


Furthermore, the u>a independence of H a p implies that the system is symmetric hy¬ 
perbolic as well as strongly hyperbolic. 

Given near Schwarzschild data for (K^, k^) on a surface S^r surrounding a Schwarzschild 
black hole, strong hyperbolicity is a sufficient condition for the system (14.471 )- (14.4811 
to produce a unique solution of the constraint problem in some neighborhood of 5?r. 
Furthermore, the problem is well-posed so that the solution depends continuously on 
the data. For linearized perturbations the solution extends globally to r = 0. 


5 Future prospects 

We have shown that the new treatment of the constraints proposed in [1] leads to 
a well-posed constraint problem for nonlinear perturbations of a Schwarzschild black 
hole in Kerr-Schild form. As is generally the case for nonlinear problems, the solution 
is only guaranteed locally in a neighborhood of the outer surface S^r on which the 
data is prescribed. The issue of a global solution to the nonlinear problem is best 
explored by numerical techniques for integrating the hyperbolic system inward along 
the p-streamlines emanating from S^r. 

The well-posedness of this problem extends to perturbations representing a Kerr 
black hole with small spin and boost. The question whether it extends further to 
a Kerr black hole with maximal spin and arbitrary boost is more complicated. Its 
resolution would depend, among other things, upon a judicious choice of the foliation 
5? p and the p-streamlines along which the evolution proceeds. This is akin to choosing 
the lapse and shift for a timelike Cauchy evolution. 

The ultimate utility of this new approach rests upon its extension to multiple 
black holes. Formally, it can be applied to the multiple black hole problem using 
a modification of the superimposed Kerr-Schild data proposed in maim, which is 
based upon the ansatz that the initial three metric for a binary black hole is given by 

hij = 5ij + 2 + 2 ( 5 . 58 ) 

where and correspond to the Kerr-Schild data for individual boosted, spin¬ 
ning black holes. In [ T6l HZ], the actual 3-metric data is only conformal to (15.5811 . 
with the conformal factor chosen to satisfy the Hamiltonian constraint. 

In our new approach to the constraints, it is possible to retain the superimposed 
Kerr-Schild initial data in its strict 4-dimensional form 

gab = gab + 2 HW + 2 hW ] , ( 5 . 59 ) 
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where £ a ^ are null with respect to the background Minkowski metric. This determines 
the initial lapse and shift as well as the initial 3-metric (j5.58j) for an evolution along the 
p-streamlines. It is is an attractive strategy because it retains much of the algebraic 
simplicity of the Kerr-Schild metric, e.g. and ‘ satisfy the background geodesic 
equation (13.9k H M and H ^ satisfy the background wave equation (13.1 IK and the 
metric can be explicitly inverted, although in a more complicated form than (I3.12K . 

Given the background metric (15.591) . the Hamiltonian constraint can be imposed 
to express the extrinsic curvature component k algebraically in terms of K \ and 
explicitly known terms via (12.71) . The extrinsic curvature components K ?J = K ij — 
1 7 ij K l i, can be freely prescribed, say, by superposition of their individual Kerr-Schild 
values. Given a suitable foliation of the initial hypersurface 5^ p and vector field p*, 
the remaining components of the extrinsic curvature data, and k, , could then 
be determined from the hyperbolic system (I2.5I) - (I2.6K obtained from the momentum 
constraint. The only data necessary are the values of K.\ and k, on a large surface 
5^r surrounding the system. The surface data for K \ and k, could be prescribed 
(again tentatively) by the superposition of their individual Kerr-Schild values. 

A major concern in such a scheme is the effect of caustics, where the ingoing 
p-streamines focus, or a crossover surface S^x where these streamlines from opposing 
points of S^r meet. For a single black hole, the streamlines can be chosen so that 
the caustics and crossovers are inside the apparent horizon, where the interior can be 
excised. However, for binary black hole data, although the caustics can be arranged 
to lie inside the black holes, the crossover surface S^x will in general span the region 
between them. In that case, unless S^x can be chosen to be a surface of reflection 
symmetry, as in the case of data for an axisymmetric head-on collision, the inward 
evolution from S^r may produce a discontinuity on S^x, he. the data induced on S^x 
may not be single-valued. 

Considerable numerical experimentation might be necessary to deal with this issue. 
The following strategy, which puts the flexibility of symmetric hyperbolic systems to 
use, is only schematic. Unlike the iterative global nature of elliptic solvers, hyperbolic 
evolution proceeds locally along the p-streamlines and can be stopped freely. This 
can be utilized to adjust the crossover surface, by numerical experimentation, so that 
it minimizes the discontinuity on S^x along each pair of intersecting p-streamlines. 
Then any discontinuity of the solution on S^x might be removed by averaging. Since 
the hyperbolic evolution of the constraint system can also proceed in the outward 
p-direction, a smooth solution, using this averaged data on can then be extended 
outward to S^r . 

The simplicity of such a scheme for binary black hole initial data is extremely 
attractive. Whether it can be successfully implemented is again a matter for numerical 
study. If such studies were indeed successful they would lead to questions of the 
utmost physical importance: Does the resulting binary black hole initial data suppress 
junk radiation? Does it give better control over the orbital and spin parameters of 
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a binary system? The sole data needed on a single large surface in the asymptotic 
region surrounding the system distinguishes this approach from other solutions to the 
constraint problem which rely on elliptic equations. Whether this feature improves 
the physical content and control of the initial data is again a matter for numerical 
investigation. 
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